The entanglement of Heisenberg chain with next-nearest-neighbor interaction 



Shi-Jian Gu^'^u Haibin Li^, You-Quan Li^, and Hai-Qing Lin^ 

^Department of Physics, The Chinese University of Hong Kong, Hong Kong, China and 
^ Zhejiang Institute of Modern Physics, Zhejiang University, Hangzhou 310027, P. R. China 

The features of the concurrences of the nearest-neighbor and the next-nearest-neighbor sites for 
one-dimensional Heisenberg model with the next-nearest-neighbor interaction are studied both at 
the ground state and finite temperatures respectively. Both concurrences are found to exhibit 
different behaviors at the ground state, which is clarified from the point of view of the correlation 
function. The threshold temperature with respective to different number of sites and the thermal 
concurrences of the system up to 12 sites are studied numerically. 
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I. INTRODUCTION 



Recently, much attention has been focused on the en- 
tanglement in the spin systems djl H S S S 

and 

indistinguishable particle systems due to the recent- 
discovered importance of entanglement in the quantum 
theory!^ 0|. These systems typically includes trans- 
verse field Ising model and anisotropic Heisenberg model. 
And the studies not only reveal the non-trial behavior of 
entanglement in the phenomenon of condensed matter 
physics, such as quantum phase transition^J ^|, but 
also shed new light on the quantum physics. However, 
most of the previous works on the spin chain mainly fo- 
cused on the model with the nearest spin exchange inter- 
action. And in most cases, the entanglement of formation 
between two spin qubits vanishes unless the two sites are 
at most next- nearest neighbors|llj. Thus it is interesting 
to investigate the problem when other kinds of interac- 
tion besides the nearest-neighbor one exist, such as next- 
nearest-neighbor interaction. This is not merely a pure 
theoretical consideration, whereas, there does exist some 
quasi-one-dimensional compounds, such as CuGeO'^|l3j| 
and NaV205^3, manifesting strong evidence of the pres- 
ence of such interaction. 

In this paper, we study the pairwise entanglement 
of the nearest-neighbor sites and of the next-nearest- 
neighbor sites in a Heisenberg chain with the next- 
nearest-neighbor exchange both at finite temperatures 
and the ground state. The entanglement of formation, 
i.e., the concurrence [Tsll is used to quantify these two 
quantities. In the following section, we first introduce 
the model, then show that the entanglement of formation 
can be calculated cither from the ground state energy at 
T = or from the partition function at finite T . In sec- 
tion ^H] we study the properties of the entanglement at 
the ground state, and discuss some special cases. Our 
results show that the presence of the interaction between 
the next-nearest-neighbor sites does not enhance the en- 



tanglement between the nearest-neighbor sites, regard- 
less it is a ferromagnetic or antiferromagnetic coupling. 
In section Hvl the thermal concurrences of a 12-site sys- 
tem as well as the threshold temperature with different 
size's systems are investigated. Finally, a brief summary 
and discussions are given in section Ivl 



II. THE MODEL FORMULATION 

The Hamiltonian of a Heisenberg chain with the next- 
nearest-neighbor interaction and periodic boundary con- 
ditions reads 

L 

H{J) = ^ [Cj (7j + i + JcTjaj+2] , 



(1) 



where L is the number of lattice sites, = {aj,a^,aj) 
denote Pauli matrices of a spin at jth site, and J is a 
dimensionless parameter characterizing the interaction 
strength between the next-nearest-neighbor sites. The 
Hamiltonian is obviously invariant under translation, and 
moreover, it has SU(2) symmetry, which manifests the 
spin conservation. Thus the reduced density matrix be- 
tween the arbitrary two sites takes the form 



Pji 




(2) 



in the standard basis |00), |01), |10), |11), and the corre- 
sponding concurrence has already been givenQ 



C = 2 max 0, |z| — Vu+u^ 



(3) 



The entities of the reduced density matrix ij^J can be 
calculated from correlation functions G, for the present 
model, they are 
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FIG. 1: The ground state correlation function G\ (left) and 
G2 (right) versus J for various sizes. The singularity at J = 
1/2 arises from the level crossing (or degeneracy). 



where G""^ 
two sites is given by 



Hence the concurrence of arbitrary 



C 



1 



0,2|G" 



G^^ - 1 



(5) 



where the SU(2) symmetry has been taken into account. 
According to statistical physics, the correlation function 
of the next-nearest-neighbor sites at finite temperatures 
is 



(6) 



where Z is the partition function and the subscript 2 
denotes distance between two sites (so does the subscript 
1 given below). At the ground state, by the Hellman- 
Feynman theorem, we have 



_ 1 dE{J) 
T=a ~ 3 dJ 



(7) 



Then the correlation function of neighboring sites is eval- 
uated as 



(8) 



where E = (H) is the internal energy of the system. 
Hence the key point is to study the two-site correlation 
function for the next-nearest neighbors 



III. GROUND STATE CONCURRENCE 

As is well known, the exact results of the Hamiltonian 
for general J have not yet been obtained except for 
some special points. At J = 0, the exact solution for the 
ground state and excited states has been well studied by 
the Bethe-ansatz method, and the correlation function 
Gf^ is simply E/3L, thus the thermal concurrence can 
be expressed in terms of the internal energy and it equals 
to 0.386 for the ground state. When J = 1/2, the ground 
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FIG. 2: The ground concurrence Ci and C2 versus J for var- 
ious sizes. Their singularities around J = 1/2 are caused by 
the level crossing (or degeneracy). 



state consists of an equal- weight superposition of the two 
nearest- neighbor valence bond statejlq: 

IVi) = [1,2][3,4]...[L-1,L] 

IV2) - [i,l][2,3] 2, L-1] (9) 

where 



[*,J]-^(|0),:|1),-|1)^|0),). (10) 

Hence the ground state concurrence can be simply writ- 
ten as 



C 



1 



2L/2 



(-1)^/2 
2L/2-2 



(11) 



which becomes 1/4 in the thermodynamic limit. 

In general case, however, we need to solve the eigen- 
value problem of the Hamiltonian for finite-size system 
numerically. It can be shown that the ground state of 
the system for J < is antiferromagnctic 17J, while for 
J > 0, many numerical results suggested that the ground 
state is antiferromagnetic for finite chain^^. Thus we 
only need to work in the invariant subspace spanned by 
those states with equal number of down spins and up 
spins, i.e., S'^.tai = 0. For example, if L = 10, the dimen- 
sion of the subspace is 252. This subspace should include 
the eigenstate with the lowest eigenvalue of the system 
due to the global SU(2) symmetry. The ground state en- 
ergy can be obtained by diagonalizing the Hamiltonian 
in this subspace. 

We show the correlation function of the nearest- 
neighbor sites Gi and that of the next-nearest-neighbor 
sites G2 m Fig. and the corresponding Gi(2) in Fig. El 
From those two figures, we see that Gi reaches a mini- 
mum at J = 0, which means that the presence of the in- 
teraction J between the next-nearest-neighbor sites does 
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FIG. 3: Threshold temperature of concurrence of the nearest- 
neighbor sites as a function of J for different size. 



not enhance the antiferromagnetic correlation between 
the nearest-neighbor sites, regardless it is a ferromag- 
netic or antiferromagnetic coupling. From this point of 
view, we deem that the presence of interaction with the 
third party generally suppresses the entanglement be- 
tween original bi-parties. This fact obviously leads to 
a maximum value of the concurrence Ci aX J — 0. More- 
over, in the region of antiferromagnetically frustrated 
coupling J > 0, the concurrence Ci is strongly affected 
by the value of J, especially around J = 1/2, as the an- 
tiferromagnetic correlation of Gi is dramatically broken 
by the frustration effect. And Ci then will be suppressed 
down to zero quickly at a critical point . On the con- 
trary, the correlation function G2 , as well as the concur- 
rence C2, between the next-nearest- neighbor sites behave 
completely in a different way. From the Hamiltonian, we 
can easily conclude that the frustrated interaction J > 
helps the formation of antiferromagnetic correlation 6*2, 
thus it is an obviously conclusion that C2 is an increas- 
ing function of J, as shown in Fig. [3 Moreover, we can 
also see from the figures that the concurrence as well as 
the correlation function are not smoothly continuous for 
all value of J. We interpret this phenomena as a conse- 
quence of the ground-state level crossing around J — \/2 
for finite-size systems. 



IV. THERMAL CONCURRENCE AND 
THRESHOLD TEMPERATURE 

At finite temperatures, the density matrix of the sys- 
tem is the summation of all states with its Boltzmann 
weight. Therefore, in the T 00 limit, the density ma- 
trix consists of an uniform distribution in the state space, 
which leads to a vanishing G^^. From Eq. ||SJ), we then 
have a zero entanglement. Thus the thermal fluctua- 
tions at high temperature always suppresses the pairwise 
entanglement, namely, the concurrence is a decreasing 
function of the temperature. So it is natural to expect 
that there exists a threshold temperature Tth at which 
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FIG. 4: Threshold temperature of concurrence between the 
next-nearest-neighbor sites as a function of J for different size. 



the concurrence vanishes. In this section, we start with 
the dependence of threshold temperature on the size of 
the system and coupling constant. 

We use numerical method to determine the threshold 
temperature for the system of L = 4, 5, . . . , 12. The re- 
sults for the concurrence of the nearest-neighbor sites are 
presented in Fig. |2| from which we observe the finite size 
effect of the threshold temperature. We find that the 
threshold temperature converge quickly as L increases, 
this is due to the fact that the concurrence here is ex- 
pressed in terms of the two-site correlation function of 
the nearest- neighbor sites, or from the other point, sim- 
ply the internal energy. And for system's energy, a small- 
ish system, such as L > 12, can well describe the physical 
properties of the thermodynamic system. So we conclude 
that the system up to L = 12 is sufficient to describe the 
threshold temperature for the infinite system. We also 
observe that the system of even (or odd) number sites 
manifest different properties. If J < 0, the threshold 
temperature of even-number-site system is larger than 
that of odd-number-site system, which is very similar to 
the results of traditional isotropic Heisenberg model [T9|. 
This is because a ferromagnetic J does not frustrated the 
tendency of singlet formation between nearest-neighbor 
sites. If J > 0, however, the situation becomes differ- 
ent, merely due to the frustration. From the left plot of 
Fig. 13 we see that the line of L = 5 crosses the other 
lines aroimd J = 0.125, so does L — 4. Consequently, 
the threshold temperature of the nearest-neighbor sites 
is not always a decreasing function of J, the exceptions 
include L — 5 and L = 7. For large system, it is expected 
to vanish around J ~ 0.7. 

In Fig^ we present our numerical result of the thresh- 
old temperature of the concurrence between the next- 
nearest-neighbor sites. We see that for the case of L = 6, 
there is no entanglement at any temperature and J. It 
is because that if L = 6, the system can be divided into 
two parts, including sites {1,3,5} and {2,4,6}, respec- 
tively. For each closed part, the frustrated interaction 
does not assist the entanglement of formation because 
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FIG. 5: Representation of thermal concurrence C between the 
nearest-neighbor sites as a function of the temperature and 
the next-nearest-neighbor interaction J. Here L = 12. 




FIG. 6: Representation of thermal concurrence C between the 
next-nearest-neighbor sites as a function of the temperature 
and the next-nearest-neighbor interaction J. Here L = 12. 



of the inside-frustration. For other cases, the thresh- 
old temperature exhibits "down-down-up-up" behavior 
as the number of sites increases, which differs from that of 
the nearest-neighbor sites completely. And it converges 
quickly as L exceeds 10. 

As we pointed out before, the entanglement of a small- 
ish system can well represent the behavior of a large 
system. We show the thermal concurrence of nearest- 



neighbor sites as a function of the temperature and the 
interaction J in Fig. [S] Generally, the thermal fluc- 
tuation suppresses the pairwise entanglement. Around 
J = 0.6, the concurrence tends to zero at the ground 
state, the thermal fluctuation may enhance it. In Fig. [SJ 
we show the results of thermal concurrence of the next- 
nearest-neighbor sites, from which we can see that there 
is no entanglement at any temperature li J < Jc- While 
at J > Jc, the thermal frustration also generally breaks 
the entanglement except in a small region of J — > Jc- 
V. SUMMARY AND DISCUSSION 

In this paper, we studied the entanglement between 
the nearest-neighbor sites as well as the next-nearest- 
neighbor sites in a Heisenberg chain with the next- 
nearest-neighbor interaction both at the ground state 
and finite temperatures. We found that the presence 
of a frustrated interaction J > will induce the entan- 
glement of formation between the next-nearest-neighbor 
sites. Thus, the entanglement in spin system stems from 
sufficiently large interaction between two sites . How- 
ever, physically, the long-range interaction usually de- 
creases as the distance increases. A typical example is 
the Haldane-Shastry model0,|^, which is modelled by 

H = Y.nmJn<ymCFm+n with J„ = Jq / 2 SIT? {niT / N) . ItS 

correlation function {aa<Tn) also decreases rapidly as n 
increases, which results in a zero concurrence even be- 
tween the nearest-neighbor sites. It is therefore believed 
that the concurrence is unwonted resources in realistic 
condensed matter. We also found that the presence of 
the interaction with third party generally suppresses the 
entanglement between original bi-parties, regardless it is 
ferromagnetic or antiferromagnetic. 

At finite temperature, we studied the entanglement by 
considering the thermal fluctuation. The threshold tem- 
perature of the entanglement was discussed for different 
system sizes. We found that Tth shows different proper- 
ties for the system with even and odd number sites, and 
it converges quickly when the number of sites exceeds 10. 
The thermal fluctuation generally suppresses the entan- 
glement at finite temperatures except in some region, i.e., 
around J = 0.6, where it may enhance the entanglement. 
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